#


import numpy as np
import matplotlib.pyplot as plt


"""
这个方程是一个关于x和y的三次方程，它的解可以通过求解以下方程得到：
x^2 + y^2 - 1 = 0 x^2 * y^3 - 1 = 0
这两个方程的解可以通过求解以下方程得到：
x^2 + y^2 - 1 = 0 y^3 = 1
第一个方程的解是x=0和y=0，第二个方程的解是y=1。
因此，这个函数的解是一个以(0, 0)为中心，以1为半径的球形区域。
"""
def heart_shape(x,y):
    return(x**2 + y**2 -1)**3 - x**2 * y** 3 <= 0

x= np.linspace(-1.5 ,1.5 , 600)
y=np.linspace(-1.5,1.5,600)
x , y=np.meshgrid(x ,y)
heart_mask = heart_shape(x,y)
plt.imshow(heart_mask, cmap='Reds', origin='lower', extent=[-1.5,1.5,-1.5, 1.5])
plt.axis('off')
plt.show()